Quantum stabilization of classically unstable plateau structures 
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Motivated by the intriguing report, in some frustrated quantum antiferromagnets, of magnetization plateaus 
whose simple collinear structure is not stabilized by an external magnetic field in the classical limit, we develop 
a semiclassical method to estimate the zero-point energy of collinear configurations even when they do not cor- 
respond to a local minimum of the classical energy. For the spin- 1/2 frustrated square-lattice antiferromagnet, 
this approach leads to the stabilization of a large 1/2 plateau with "up-up-up-down" structure for J2/J1 > 1/2, 
in agreement with exact diagonalization results, while for the spin- 1/2 anisotropic triangular antiferromagnet, 
it predicts that the 1/3 plateau with "up-up-down" structure is stable far from the isotropic point, in agreement 
with the properties of Cs 2 CuBr 4 . 



Introduction. — Frustration is responsible for the emergence 
of several remarkable properties in quantum magnets, rang- 
ing from rather exotic types of order such as quadrupolar or 
nematic order to resonating valence bond or algebraic spin 
liquids [1]. In the presence of an external field, frustration is 
also known to be at the origin of several types of accidents in 
the magnetization curve, including kinks, jumps and plateaus. 
Of all these remarkable features, magnetization plateaus at ra- 
tional value of the magnetization are probably the best docu- 
mented ones experimentally, and their theory is likewise quite 
advanced. Following the terminology of Hida and Affleck 
[2], two kinds of plateaus have been identified [3]: 'classical' 
plateaus [4-6], whose structure has a simple classical analog 
with spins up or down along the external field, and 'quantum' 
plateaus [7-11], which have no classical analog and corre- 
spond to a Wigner crystal of triplets in a sea of singlets. In the 
case of quantum plateaus, the mechanism is clear: frustration 
reduces the kinetic energy of triplets, resulting in a crystal- 
lization at commensurate densities. The main open problem 
is to be predictive for high commensurability plateaus since 
it requires a precise knowledge of the long-range part of the 
triplet-triplet interaction. 

By contrast, and somehow surprisingly, the theory of clas- 
sical plateaus is not complete yet. The paradigmatic example 
of a classical plateau is the 1/3 magnetization plateau of the 
Heisenberg antiferromagnet on a triangular lattice, studied by 
Chubukov and Golosov [5] in the context of a l/S expan- 
sion. In this system the three sublattice up-up-down (uud) 
structure appears classically atH- H sal /3, and quantum fluc- 
tuations stabilize this uud state in a finite field range around 
i/ sa t/3, leading to the 1/3 plateau. The basic qualitative idea in 
the spirit of the order by disorder is that collinear configura- 
tions often have a softer spectrum and, hence, a smaller zero- 
point energy [12, 13]. The prediction of the 1/3 plateau has 
been confirmed by exact diagonalization of finite clusters for 
5=1/2 and 1 [14], and the theory of Chubukov and Golosov 
can be extended to all cases where a collinear state is classi- 
cally stabilized for a certain field. 

There are cases, however, where a classical plateau has 
been suggested to exist although the collinear structure sta- 



bilized for quantum spins is not the ground state for classical 
spins, the classical ground state in the appropriate field range 
being in general a non-coplanar structure. This is for instance 
the case of the spin- 1/2 J\-J2 model on the square lattice, for 
which exact diagonalizations have revealed the presence of a 
four sublattice up-up-up-down (uuud) 1/2 plateau in a param- 
eter range where the classical ground state has a canted stripe 
structure (see below). Another example is the 1/3 plateau of 
the Heisenberg model on the anisotropic triangular lattice, a 
model relevant to the compound CS2CUB14. To develop a gen- 
eral theory of classical plateaus in that situation remains the 
main open issue in the field. 

The goal of this Letter is to develop such a theory. For 
that purpose, we start with a general Heisenberg model in an 
external field defined by the spin Hamiltonian 

'H = Y J J ij S i Sj-HY J S^ (1) 

<ij> i 

and show how to estimate the zero-point energy of collinear 
states even if they do not minimize the classical energy. More 
precisely, we derive an upper bound of this energy to order 
1/5. If the energy of a collinear state estimated in this way 
is lower than that of the classical ground state (including the 
zero-point energy), then the collinear state must be the ground 
state since the energy used for the comparison is an upper 
bound. Correspondingly, the quantum antiferromagnet ex- 
hibits a magnetization plateau in a certain field range, which 
may generally exceed our conservative theoretical estimate 
based on the upper energy-bound. We apply this approach 
to the J\-J2 model on a square lattice and to the Heisenberg 
antiferomagnet on the anisotropic triangular lattice, with re- 
sults in remarkable agreement with existing numerical data 
forS = l/2. 

General formalism. — We propose a simple method to test 
for the stability of plateau structures in the framework of 
the linear spin-wave theory. First of all, we note that, since 
plateau structures are collinear, the fluctuation Hamiltonian 
around such structures in terms of Holstein-Primakoff bosons 
does not contain linear bosonic terms, and that the first rele- 
vant terms in its l/S expansion are quadratic. Since plateau 
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structures are classical minima of the energy only at specific 
values of couplings and magnetic field (if any) the corre- 
sponding harmonic fluctuation Hamiltonians are positive def- 
inite only at these specific points. Away from such points the 
correction to the classical energy cannot be straightforwardly 
computed. The fact that the spectrum is not well defined stems 
from the harmonic approximation. If the plateau state is to 
become the true quantum ground state, higher order terms in 
the spin-wave expansion must produce an excitation spectrum 
with positive frequencies. This approach requires to calcu- 
late nonlinear quantum corrections to the spectrum as done in 
Ref. [5] and in more recent studies [15, 16]. This is a rather 
cumbersome procedure, and it would be useful to have a sim- 
pler approach to determine if there is a plateau, and to esti- 
mate its width to the lowest-order in l/S . Besides, the calcu- 
lation of the excitation spectrum for the plateau region allows 
to identify only second-order transitions, while in many ex- 
perimental and model examples the transitions at the plateau 
boundaries are of the first-order and, therefore, require a full 
energetic comparison. 

The proposed method to obtain a well-defined spectrum 
around a state which is not a classical ground state consists 
in adding a staggered-field term 

V = 5}](S-S?) (2) 

i 

to the harmonic Hamiltonian. A similar approach has been 
introduced in a different context in [17, 18]. On each lattice 
site i, the staggered field 5 > is oriented in the direction z, of 
the corresponding classical spin. The extra term (2) amounts 
to a shift of the chemical potential of the Holstein-Primakoff 
bosons V = 5Yji a]ai and yields a positive contribution to the 
spin-wave Hamiltonian. The magnitude of 5 is adjusted to en- 
sure that the harmonic Hamiltonian is positive definite. The 
resulting spectrum obtained with the help of the Bogolyubov 
transformation has real and positive frequencies. The advan- 
tages of this variational approach are the following: the ad- 
dition of V to the Hamiltonian does not change the classical 
energy of the trial state and allows to obtain physically mean- 
ingful dispersion relations. Furthermore, since the expectation 
value of V is strictly positive, the computed energy correction 
provides an upper bound for the energy of the plateau state. It 
should be noted that the suggested method is valid not only for 
collinear spin structures but can be extended to all structures 
that are saddle points of the classical energy, see [17, 18] and 
below. 

J\-J2 model. — For the frustrated square-lattice antiferro- 
magnet the exchange interaction constants are 7 (/ = J\ and 
J2 for nearest and second-nearest neighbors, respectively. In 
zero field the classical ground state of the model is a helix with 
the ordering wave vector given by the minimum of the Fourier 
transform of the coupling interaction J q = 4/iy q + 4/2^ 
with y q = (cos q x + cos q y )/2 and rj q = cosq x cosq y . For 
J2/J1 < 1/2, the minimum corresponds to qjy = (n, tt), i.e. 
to Neel order. In the opposite case J2/J1 > 1/2, the order by 
disorder mechanism selects collinear striped structures with 



ordering wave vectors q$ = (n, 0) or (0, n) [13]. The point 
J2lJ\ - 1/2 is highly degenerate since 7 q is minimal along 
the lines q x = n and q y = n. In the presence of a magnetic 
field both the Neel and the stripe structure are canted with 
a uniform spin component in the field direction. The cant- 
ing angle 8 measured with respect to the z axis is given by 
cos On = H/8J1S and cos# s = ///(47, + 8J 2 )S for the two 
states. 

The analysis of classical spin configurations in a magnetic 
field suggests the appearance of a 1/2-magnetization plateau 
with a four-sublattice uuud structure for the strongly frustrated 
point J2/J\ =1/2 [6]. The conclusion has been supported by 
exact diagonalizations of finite clusters, though numerically 
the plateau extends well into the classically unstable region 
hlh > I/ 2 with the largest width at J 2 /J\ ~ 0.6. More- 
over, the linear spin-wave calculation for J2/J1 = 1/2 shows 
that for this ratio of coupling constants the canted Neel state 
wins over the collinear plateau state [19] leaving an apparent 
problem with reconciling numerical and analytical results. 

We now investigate the appearance of the 1/2- 
magnetization plateau for the J1-J2 model using the 
variational harmonic theory outlined above. In the semiclas- 
sical approach deviations from the classical configuration are 
expressed as Holstein Primakoff bosons [20]. In the harmonic 
approximation the bosonic Hamiltonian can be split into three 
contributions: 

<H = <H (SS) + ( H {1) +"K (2) , (3) 

where '7Y (0) is the classical energy of the system. and 
'MS® respectively contain only terms which are linear and 
quadratic in boson operators. The uuud state, being a collinear 
state, is such that no linear terms appear when performing the 
HP transformation, hence f HS v > = for this state. In non 
collinear structures r H m is proportional to the derivative of 
the classical energy with respect to spin orientations. It van- 
ishes for the canted Neel and canted stripe structures over the 
entire parameter range since both are extrema of the classical 
energy. In fact the canted Neel (stripe) structure is a saddle 
point of the classical energy for H < i/ sat and J2IJ1 > 1/2 
(respectively J2IJ1 < 1/2). 

The general structure of the quadratic bosonic Hamiltonian 
that describes harmonic fluctuations around these various con- 
figurations is given by 

ti = NE C \ + 1 ~Yj KM k a k - A k ] , (4) 
1 k 

where E c \ is the classical energy per site of the state around 
which fluctuations are considered. For the canted Neel and 
the canted stripe states a k = (a k ,fl- k ) and M k is the 2x2 
matrix 

Coefficients A k (q) and B k (q) for Neel and striped structures 
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FIG. 1 : Semiclassical phase diagram for the spin 1/2 Ji — Ji Heisen- 
berg model in magnetic field. The uuud structure is stabilized by 
fluctuations over a wide parameter range. Dashed lines correspond 
to the phase diagram for 5 = 1 and dotted lines to the classical phase 
boundaries. The shaded area represents schematically the gapped 
singlet phase for the spin- 1/2 model. 

are listed below: 

A k (q N ) = 4J l S(l+y k cos 2 N )-4J 2 S(l-T ]k ), 
B k (q N ) = ~4JiSy k sin 2 N (6) 

and 

^k(Qs) = 4-J2S(l+ Tj k cos 2 8s) + 2J\S(cosk y + cos 2 6s cosk^, 
B k (q s ) = -25sin 2 s (y 1 cosfc A . + 27 2 77 k ) (7) 

where we have used q$ = (n, 0). The additional constants 
A^ and A£ in Eq. (4) are given respectively by Ak(q^) and 
A k (q s )- 

The uuud state has a four-site unit cell, and a k denotes 
(a[ k , ■ ■ • ,a\^,a\- k , ■ ■ ■ ,fl4,-k) with M k being the 8x8 matrix 
obtained from (5) by substituting 

' A k E k G* \ ( -f k \ 

E* B k F k OH* 

Ak ~ C k ' ° k ~ -F* H k -E k 
\G k F* A k ) { -E* J 

(8) 

with coefficients 

A k = -4J 2 S + H , B k = -4(/i - J 2 )S + H , 

C k = 4(7! + J 2 )S -H, E k = JiSt^ , 

F k = JiSTk x , G k - J 2 S T-k x T-k } , H k = -J 2 S T- kx T ky , 

(9) 

where T k = (1 + e~ ak ). The constant Aj[" iK/ is given by 
2A k + B k + C k . When the quadratic form of Eq. (4) is positive 
definite, it can be diagonalized by the standard Bogolyubov 
transformation allowing to compute the quantum corrections 
to the classical energy. For parameters for which the matrix 
M k is not positive definite we add to the Hamiltonian the term 
V defined in Eq. (2). From the expression of V in terms of HP 



bosons it is clear that it leaves < K (0) and 1ri m unchanged while 
its effect on the bosonic Hamiltonian is to increase Ak and all 
diagonal elements of M k by 6/2. The field 6 is adjusted to the 
minimal value which is sufficient to make Mk positive definite 
over the entire Brillouin zone. 

The phase diagram obtained by comparing the ground-state 
energies for three relevant spin structures is presented in Fig. 1 
for 5 = 1/2 and 5 = 1. The 1/2-magnetization plateau is sta- 
bilized by quantum fluctuations over a wide range of param- 
eters deep into the classically forbidden region J 2 /J\ > 1/2, 
though it remains energetically unfavorable at J 2 /Ji = 1/2. 
The width and position of the plateau are in good agree- 
ment with the exact diagonalization results of finite clusters 
with up to = 36 sites [6], Figure 2 shows the magnetiza- 
tion curves for several ratios J 2 IJ\. The magnetization curve 
for J 2 /Ji = 0.6 with a large magnetization jump below the 
plateau and a much smaller anomaly above the plateau is in 
good correspondence with the numerical data for the same 
coupling ratio [6]. For J 2 /J\ close to 1/2, there is in addi- 
tion a competition between the canted Neel and the canted 
stripe states. The Neel state has a softer spectrum than the 
stripe state and is stabilized beyond its classical boundary. 
This leads to an additional transition from the canted Neel 
state into the canted stripe structure which shows up as a small 
jump either above (J 2 /Ji = 0.525) or below (J 2 /J\ = 0.55) 
the 1/2-plateau. 

It should be pointed out that other states than those consid- 
ered may be stabilized. One possible candidate, at the upper 
edge of the plateau, is the coplanar four sublattice state having 
three classical spins parallel and the remaining spin pointing 
in a different direction. This state can be naturally connected 
to the uuud state and is the analog of the state stabilized above 
the plateau in the isotropic triangular lattice. However, such 
a structure could not be investigated in our linear spin-wave 
approach since it is not a saddle point of the classical energy. 
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FIG. 2: Magnetization curves of the spin 1/2 J\ - J 2 model for dif- 
ferent ratios Ji/Ji obtained in the variational spin wave approach. 

Anisotropic triangular lattice. — We now consider a sec- 
ond example of classically unstable magnetization plateau, 
the nearest-neighbor Heisenberg antiferromagnet on an or- 
thorhombically distorted (anisotropic) triangular lattice. In 
this model spins are coupled by J/j = J along horizontal 
chains and by 7 (/ = J' on zig-zag interchain bonds, see 
Fig 3(a). The spin-1/2 model is relevant for Cs2CuBr4 [21- 
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23], which has a 1/3 magnetization plateau although, with 
J' ~ 0.757, it is quite far form the isotropic limit. The robust- 
ness of the uud plateau in the J-J' model has been studied nu- 
merically [24] and analytically [15]. Nevertheless the extent 
of the plateau state around H/H sat ~ 1/3 and J' /J ~ 1 and the 
nature of the states adjacent to the plateau region are still open 
questions. Furthermore, the stability method employed by Al- 
icea et al. [15] allows to identify only second-order transitions 
out of the plateau state, while experiments typically find first- 
order transitions [23]. We shall see below that this fact finds a 
natural explanation in our theoretical approach. 

The Fourier transform of the coupling interaction in the 
triangular-lattice, see Fig. 3(a), is given by J q - 2[/cosqa + 
/' cos qb + /' cos q(a - b)]. In zero field the classical ground 
state is a helical spin structure whose ordering wavevector Q 
minimizes J q . In the isotropic case J-J' this yields the well- 
known 120°spin structure. In the presence of a magnetic field 
the classical energy is minimized for canted helices or um- 
brella configurations, see Fig. 3(c), which have helical order 
in the xy plane and uniform spin component in the field direc- 
tion. The canting angle of the helical structure measured with 
respect to the z axis is given by cos 6h = H/(Jq-Jq)S . For the 
isotropic point J' /J = 1 the canted helical state is degenerate 
with the coplanar Y- and V-type structures, see Fig. 3(b). Ex- 
isting linear spin-wave calculations indicate that the coplanar 
structure is selected over the non coplanar one in the isotropic 
lattice and that the uud structure, classically stable at the field 
i/ sa t/3, is stabilized by fluctuations over a finite field range [5]. 

In the following we address the problem of the plateau sta- 
bility for the anisotropic triangular lattice model by compar- 
ing ground-state energies for the canted helical state, the uud 
structure and the two 3-sublattice planar states as a function of 
J' I J and magnetic field. Away from J = J', the 3-sublattice 
planar structures turn out to be, like the iW-state away from 
i/ sa t/3 and J = J', saddle points but not local minima of the 
classical energy, and to compute their zero-point energies, we 
use the variational spin-wave approach suggested above and 
add a staggered local field to the Hamiltonian [29]. 

The resulting phase diagram is shown in Fig. 4. The uud 
plateau state is stabilized well beyond the isotropic limit and 
extends over the range 0.5 <J' /J < 1.5. Coplanar states are 
stabilized above and below the magnetization plateau with the 
exception of the plateau edges, where we find direct first-order 
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FIG. 3: a) Anisotropic triangular lattice and basis used in our calcula- 
tions, b) 3- sublattice planar structure as a function of the magnetic 
field, c) Example of non coplanar canted helix. 
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FIG. 4: Phase diagram of the spin- 1/2 anisotropic triangular lattice 
in magnetic field. Y and V regions denote 3-sublattice planar states. 
The dashed line is the classical saturation field. The gray shading 
denotes regions where other phases than the canted helical states may 
be expected. 



transitions from the uud state into the canted helical structure. 
From the energetic comparison it appears that the uud state 
does not extend into the Y-state region, so that the correspond- 
ing portion of the lower boundary of the plateau is perfectly 
linear (see Fig. 4). This is almost certainly an artefact of the 
method, which only gives an upper bound to the energy of the 
plateau, and the extent of the plateau is probably significantly 
larger. In fact, at J — J', we obtain a plateau width which is 
only half that predicted in Refs. [5, 14]. We also note that the 
only coplanar states considered in our calculation are Y- and 
V-type structures, while for a substantial mismatch between J 
and incommensurate coplanar structures may be also stabi- 
lized by quantum fluctuations. The variational spin-wave ap- 
proach is not well-suited for treating them and we only remark 
that they may appear on the phase diagram at the expense of 
the canted helical structure. Finally, the gray shading in Fig. 4 
indicates regions, where new quantum phases are expected. In 
fact, in zero field, theoretical and numerical approaches point 
to collinear spin correlations for weakly coupled chains, [25- 
27], while in the limit of strong interchain couplings the AFM 
Neel state should be stable down to J'/J « 1.5 [27]. 

As compared to those of Alicea et al, who also predicted 
an extended plateau region for small distortions [(1 - J'/J) 2 < 
0.3] for the 5=1/2 case [15], our results bring in a num- 
ber of new insights. In the first place, the symmetry between 
J'/J < 1 and J'/J > 1 is lost. Second, a transition out of 
the plateau into the canted helical states is clearly present. 
Finally, for J'/J = 0.75 relevant for Cs 2 CuBr 4 [21-23], we 
find a magnetization plateau width AH/H yit « 0.106, signifi- 
cantly larger than the experimental value AH/H sdt ss 0.052 in 
Cs2CuBr4. Since in most cases our approach underestimates 
the plateau width, the difference must be due to additional 
effects not included in the anisotropic model, for instance 
the competition between quantum effects and Dzyaloshinskii- 
Moriya interactions [28]. 

Conclusion. — We have developed a general method to in- 
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vestigate the stabilization of classical magnetization plateaus 
in cases where the corresponding configuration is not a mini- 
mum of the classical energy. This method is extremely simple 
since it only relies on the diagonalization of quadratic bosonic 
Hamiltonians and does not require to go beyond linear spin- 
wave theory, yet it appears to give remarkably accurate results, 
even for spin 1 /2. This has been demonstrated in two cases 
of current interest, the /1-/2 Heisenberg model on the square 
lattice and the Heisenberg model on the anisotropic triangu- 
lar lattice, for which it predicts that plateaus at magnetization 
1/2 and 1/3 respectively are stabilized over a wide range of 
parameters. 
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MaNEP. 
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SUPPLEMENTAL MATERIAL: LINEAR SPIN WAVE 
THEORY FOR THE ANISOTROPIC TRIANGULAR 
LATTICE AFM IN A MAGNETIC FIELD 

The general structure of the bosonic fluctuation Hamilto- 
nian around a given classical state is given in Eq. (4) of the 
main text. In this supplemental material we present the spe- 
cific expressions for the canted helical state, the uud state and 
the 3- sublattice coplanar structures. 

Canted helical state — The classical energy per site of 
the canted helical state expressed as a function of the canting 
angle is given by 

E% = S 2 (J cos 2 8 + Jq sin 2 0)/2 - HS cos 6. (10) 

It is minimal for cos0 H = H/(Jq - Jq)S . The term a k in 
Eq. (4) of the main text denotes a_ k ) and M k is the 2 x 2 
matrix 



M k (Q) 



A k (Q) + C k (Q) B k (Q) 
fik(Q) A k (Q) - C k (Q) 



(11) 



with coefficients 



A k (Q) = S (-J Q + i(cos 2 + l)(/ k+Q + / k _ Q ) + \ sin 2 0J k ) , 
Sk(Q) = S (i(cos 2 - l)(7 k+ Q + 4-q) + \ sin 2 0J k ) , 
C k (Q) = 5cos0(7 k+Q -/ k _ Q )/2. 

(12) 

The additional term A k in Eq. (4) of the main text is given by 
A k = -S/ Q . 

3-sublattice coplanar states — The classical energy per 
site of any 3- sublattice coplanar structure is given by 

^coplamu _ g2(y + 2/')(cOS 0^2 + COS 0^3 + COS »2,3)/3 

-SH(cos a\ + cos + cos ff3)/3 

(13) 

where atj = cti-aj are the spin orientations measured with re- 
spect to the field direction. The angles a, minimizing Eq. (13) 
are given by 



aY = 7T, cos aY = - ( — + 1 1 , at = -a! (14) 

1 2 2Usat/3 ' 
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for < H < ff sat /3 and 



where yP^ and y'f' are given by 



(2) 



H 



cos a. 



2H^ y H 2 

H L H 2 



cos as 



■ cos as 



4H s; 



3 + 



H 2 



(15) 



for fields in the range H sat /3 < H < Z/ sat , where H sat is the 
saturation field i/ sat = 3(7 + 2J')S . 

Since the states considered have three sites per unit 
cell the term a k in Eq. (4) of the main text denotes 

( a I l ' a I 2' a I 3' a -k,i ' fl -k,2> fl -k,3) an d Mk is the 6x6 matrix 



M. 



Y,V 



A D* H* E* /* 

D k B F* E k G* 

#k fk C Ik Gr 

£* /* A D* H* 

E k G* D k B F* 

Ik Gk Hk Fk C , 



(16) 



^(2) _ y g ik(-a+2b) + y/ g /k(a+b) A + g /k(-a+2b)^ _ 



(18) 



with 



A = [-S(J + 2/')( cos ffi 2 + cos a\ t x) + Hcos ct[] , 

B = \—S(J + 2/')( cos (X\,2 + COS a2,3) + HCOS O^] , 
C = [S(J + 2y')( cos <*2,3 + COS a l,3) + ^ cos a 3] > 

D k = S(cos aia + l)r^/2, £ k = 5(cosa li2 - 1)7^/2, 
F k = S (cos a 2 ,3 + 1)7^/2, G k = 5 (cos a 2 ,3 - 1)7^/2, 
H k = 5(cosa lt3 + l)r k 2) /2, 4 = S(cosari >3 - l)7 k 2) /2. 



(17) 



The additional term Ak in Eq. (4) of the main text is given by 
A k = A + B + C. 



uud state — The uud state also belongs to the family of 
coplanar states. Hence the expressions for the classical en- 
ergy and for the coefficients of the fluctuation Hamiltonian 
(Eqs. (13), (16) and (17)) obtained previously can be applied 
to the uud state if one replaces a\ = n and ar 2 — — 0. 



